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Abstract 

> 

\^ ' An effective cylindrical cloak may be conceptualized as an assembly of adjacent local neighbour- 

fT^ , hoods, each of which is made from a homogenized composite material (HCM). The HCM is required to 

0^ ' be a certain uniaxial dielectric-magnetic material, characterized by positive-definite constitutive dyadics. 

in ■ It can arise from the homogenization of remarkably simple component materials, such as two isotropic 

QQ . dielectric-magnetic materials, randomly distributed as oriented spheroidal particles. By carefully control- 

' ling the spheroidal shape of the component particles, a high degree of HCM anisotropy may be achieved, 

which is necessary for the cloaking effect to be realized. The inverse Bruggeman formalism can provide 
estimates of the shape and constitutive parameters for the component materials, as well as their volume 
^ fractions. 
• ^ - 

X 

' Keywords: cloak, geometric optics, Bruggeman formalism, inverse homogenization 

1 Introduction 

The problem of light propagation in a homogeneous material occupying a region of complex shape can 
be reversibly transformed into the problem of light propagation in an inhomogeneous material occupying 
a region of simple shape. This observation has been appreciated since Tamm's pioneering work on the 
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electromagnetics of curved spacetime in the 1920s [1, 2], and it also provides the foundation for the well- 
established C method in grating analysis, developed by Chandezon and colleagues [3, 4] and applied in 
the context of negative refraction [5]. Additionally, it underpins recent proposals of electromagnetic cloaks, 
which supposedly render objects contained inside the cloaks almost invisible to an external observer [6] . These 
cloaks are made of composite materials with judiciously-engineered microstructures or nanostructures. 

A variety of theoretical schemes for achieving electromagnetic cloaking have been put forward [7, 8, 9, 
10]. Furthermore, partial cloaking in the microwave frequency range — at a single frequency for a single 
polarization — was observed in the laboratory five years ago for a cloak made from an assembly of split-ring 
resonators [11]. More recently, there have been further experimental reports of cloaking involving tapered 
waveguides [12] and plasmonic composite materials [13], for examples. 

Designers of composite materials face severe challenges when attempting to realize a cloak. Typically, 
a high degree of anisotropy and a high degree of inhomogeneity are needed. In order to cope with these 
requirements, cloak designs are generally based on rather complex microstructures or nanostructures. How- 
ever, some simplifications have emerged very recently: A lamination process can be implemented in order to 
realize a wide range of dielectric and magnetic constitutive parameter values, such as may be required for 
an electromagnetic cloak [14]. Also, mention should be made of a relatively simple cloak design based on an 
array of elliptical rods proposed by Gao et al. [15]. 

The aim of realizing a cloak using a simple microstructurc or nanostructurc provides the motivation 
for our study. In the following we show how the homogenization of a random distribution of oriented 
spheroidal particles, made from isotropic dielectric-magnetic materials, can give rise to a well-functioning 
quasi-monochromatic cloak, as demonstrated by a study of ray trajectories. In addition, we show how the 
constitutive and morphological parameters of the component materials, as well as their volume fractions, 
can be estimated using a process of inverse homogenization. 

In the following, vectors are underlined, with the ' symbol signifying a unit vector. With respect to the 
Cartesian unit basis vectors {x, the position vector is expressed asr = xx + yy + zz_. Dyadics are 
double underlined. The permittivity and permeability of free space in the absence of a gravitational field 
are denoted as Eq and /Ho, respectively. 



2 Cloak inspired by a cosmic string 

The object to be cloaked is contained within an empty cylinder which is aligned with the z axis and has 
radius p — pc > where the radial distance p ~ \J + y^- The constitutive relations of the inhomogeneous 
material from which the cloak is composed are 



Rir) = e„2ip) ■ E{r) 

where the diagonal 3x3 dyadic 



P>Pc, (1) 



= 711 + 11) + 722 (2) 

has components 



711 (P) 



722 (P) 



p 
p 



}■ (3) 



p-1 

The constitutive parameters 711 and 722 arc plotted versus radial distance p in Fig. 1. Both constitutive 
dyadics are obviously positive definite for all p > 1. In the limit p ^ 1, we see that 711 becomes null- valued 
while 722 becomes unbounded. On the other hand, both 711 and 722 rapidly approach unity as p — ?> 00. 

In fact, the structure of our proposed cloak is identical to the Tamm medium which represents the 
exterior region of a (non-spinning) cosmic string [16]. In a geometric-optics study of this Tamm medium, 
it was demonstrated that the interior of the cosmic string is almost entirely invisible to a distant observer 
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[17]. Furthermore, the similarity between the constitutive relations (l)-(3) and those of the experimentally- 
realized cloak of Schurig et al. [11] is noted. 



3 Ray trajectories 

A study of light rays passing through the cloak specified by the constitutive relations (l)-(3) was performed, 
within the geometric optics regime. As our ray-tracing methodology is described in detail elsewhere [17, 18], 
only a brief outline is provided here. For a quasi-planewave with relative wavevector k, the quantity 



represents a suitable Hamiltonian function. With the paramcterizations rir) s-nd ^(t), ray trajectories are 
governed by the coupled vector differential equations^ [19, 20] 



with their direction being given by Vi/H, which is aligned with the time-averaged Poynting vector. Standard 
numerical methods, such as the Runge-Kutta method, can be used to solve the system (5) for r(T) and fc(T), 
once appropriate initial conditions r(0) and fc(0) have been specified. 

An example of planar ray trajectories is illustrated in Fig. 2 for pc = 1.1. Yellow shading indicates the 
disc p = 1. The rays start at equally-spaced locations in the first quadrant of xy plane along a line parallel 
to y = —X, with fc(0) directed along (—1, —1, 0). The ray trajectories remain in the xy plane. The rays skirt 
around the outside of the circle p = Pc, such that within a few radiuses from the centre of the circle there is 
minimal deflection of the rays from their original straight line trajectories. Thus the contents of the cloak 
are almost entirely invisible to a distant observer located in the third quadrant of the xy plane. 

An example of ray trajectories in three dimensions is represented in Fig. 3. Here rays start at equally- 
spaced locations within a plane parallel to y = —x in the octant {x > 0,y > 0, z > 0}, with fc(0) directed 
along (—1, —1, —1). The rays do not remain in a single plane in this case, but the deviation from the plane 
is relatively small, especially at distances beyond a few radiuses from the cloak's axis. 

By considering values of pc smaller than 1.1, the degree of cloaking is improved very little as perceived 
by an observer more than a few radiuses away from the cloak's axis. However, by increasing the value of pc 
from 1.1 the degree of cloaking quite rapidly diminishes. 

4 Cloak as a homogenized composite material 

Let us now turn to the issue of realizing the cloak specified by the constitutive relations (l)-(3). It was 
recently demonstrated that certain uniaxial dielectric-magnetic Tamm mediums can be envisaged as homog- 
enized composite materials (HCMs) that can be fabricated from relatively simple component materials [21]. 
A similar approach is taken here. 

Before turning to the homogenization process itself, we must first deal with the inhomogeneity of the 
proposed cloak. This may be catered for by subdividing the region of space occupied by the cloak into 
local neighbourhoods which are sufficiently small as to be considered approximately homogeneous. Then 
the inverse homogenization procedure described later in this section can be applied locally. This piecewise 
homogeneous approach is documented in detail elsewhere for general Tamm mediums [22]. 

Now we turn to the homogenization of simple-component composite materials, with our aim being to spec- 
ify a HCM whose constitutive parameters coincide with the proposed cloak. As our theoretical foundation, 
we rely on the well-established Bruggcman formalism [23, 24]. 

*Here the shorthand Vq = xd/dqx + y djdqy + z_d/dqz for q = x + qy y + i is adopted. 



Ti = det 7(p) - k ■ j{p) ■ k 



(4) 




(5) 
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Suppose we consider the homogenization of a composite material comprising two isotropic component 
materials, one being a dielectric-magnetic material with relative permittivity Ca and relative permeability 
/^a and the other likewise being a dielectric-magnetic material but with relative permittivity ef, and relative 
permeability /^b- The component materials a and h are randomly distributed with respective volume fractions 
fa and fb = i—fa- Each component material is distributed as spheroidal particles. The axis of these spheroids 
for both component materials is assumed to be aligned with the symmetry axis of 7, namely the y axis. 
Thus, the vector 

r_g = V( {x_x_ + Uii/y + z z) • r , {i = a, b) , (6) 

prescribes the surface of each spheroid relative to its centre. Here the vector r_ prescribes the surface of the 
unit sphere, while the parameter Vi > is el linear measure of size which is taken to be much smaller than 
the electromagnetic wavelengths involved. Prolate spheroids are characterized hy U( > 1 whereas oblate 
spheroids are characterized by Ui < 1. 

In accordance with the Bruggeman formalism, the corresponding HCM is a uniaxial dielectric-magnetic 
material. Its relative permittivity dyadic e^'' and relative permeability dyadic fi^^' have the form 



'11 



ixx + zz)+T.22''yy, (r = e,/i). (7) 



For the case of uniaxial dielectric-magnetic HCMs such as is involved here, full details of the numerical 
process of computing the dyadics e^^ and /i^*" are available elsewhere [21, 23]. 

Usually, homogenization formalisms are implemented in the forward sense, whereby the constitutive 
parameters of the HCM are estimated from a knowledge of the constitutive and morphological parameters of 
their component materials. However, our goal here is to estimate constitutive and morphological parameters 
of the component materials which give rise to a HCM such that = /_tf/" =711 and = 1122 = 722- 
In order to do so, an inverse implementation of the Bruggeman formalism is required. We note that formal 
expressions of the inverse Bruggeman formalism have been established [25], but in some instances these 
formal expressions can be ill-defined [26]. In practice, direct numerical methods may be more effective 
in implementing the inverse formalism [27]. A note of caution should be added concerning the inverse 
Bruggeman formalism: certain constitutive parameter regimes have been identified as problematic [28] , but 
these regimes do not overlap with those considered here. 

As a representative example, let us focus on a particular inverse homogenization scenario whilst noting 
that other scenarios are also feasible. For convenience, we seek component materials which are composed 
of spheroidal particles all of the same shape; i.e., Ua = Ub- Furthermore, we shall assume that the relative 
permittivity and relative permeability parameters of the component materials are the same; i.e., Ca = Sb and 
fXa = fJ-b- This coincidence of constitutive parameters can be straightforwardly engineered if the component 
materials themselves are envisaged as HCMs, as has been demonstrated previously [21]. Now, assuming 
that the relative permittivity and relative permeability of component material a are fixed, we implement 
the inverse Bruggeman formalism in order to estimate the volume faction, shape parameter and constitutive 
parameters for component material b. 

The inversion of the Bruggeman formalism may carried out as follows. Let {efi^, £^2*^, Mn^, ^^2^} denote 
the forward Bruggeman estimates of the HCM's relative permittivity and relative permeability parameters, 
computed for physically-plausible ranges of the parameters fb, Ub and eb, namely fb G {fb,fb^), Ub € 
{U-,U+) and eb G (efc-,e^). Then: 

(i) Fix Ub = {U- + U+) /2 and €b = (e" + e+) /2. For aU values of fb G {fbJb)^ ^nd the value fl for 
which the quantity 



\ / ^Br \ / " Br \ / " Br 

711 \ , / £22 - 722 \ , / Mil - 711 \ , / M22 - 722 



711 / \ 722 / V 711 / V 722 



1/2 



(8) 



is minimized. 
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(ii) Fix fb = fl and ej, = (e^ + e^) /2. For all values of Ub E (U^^ , U^) , find the value for which A is 
minimized. 

(iii) Fix fb = fl and Ub = U^. For all values of et e (ej^, e'^), find the value e]^ for which A is minimized. 

The steps (i)-(iii) are repeated, using and as the fixed values of Ub and in step (i), and el as the 
fixed value of eb in step (ii), until A becomes sufficiently small. 

Estimates of fb, Ub and eb computed using this scheme are presented in Fig. 4 for the case pc = 1.1. 
These quantities are plotted versus the radial distance p. The relative permittivities chosen for component 
material a were: ea = 0.01 for the range Pc < P < 1-5; Ca = 0.1 for the range 1.5 < p < 2.5; and ea = 0.3 for 
the range 2.5 < p < 5. For all the numerical results presented in Figs. 4, the degree of convergence of the 
numerical schemes was < 0.1%. While the constitutive parameters chosen for component material a are not 
commonplace, reports of materials with relative permittivities and relative permeabilities close to zero are 
to be found in the recent literature [29, 30, 31]. Materials with extreme properties are now being considered 
as well [32]. Furthermore, the constitutive parameters for component material b represented in Fig. 4 are 
not at all exotic. 

5 Closing remarks 

We have shown theoretically that the homogenization of a composite material comprising two simple isotropic 
component materials can result in an HCM which can be used to create a well-functioning cylindrical cloak. 
By carefully controlling the spheroidal shape of the particles which make up the two component materials, 
the high degree of HCM anisotropy which is essential for cloaking can be achieved. The inverse Bruggeman 
formalism provides estimates of the shapes needed for the component particles, as well as their appropriate 
constitutive parameters and volume fractions. Provided the component materials are weakly dispersive in 
their constitutive parameters, the proposed cloak can function in a spectral regime of appreciable bandwidth. 

The homogenization scenario presented here is by no means the only way of realizing a uniaxial dielectric- 
magnetic HCM of the required form. For example, the two component materials could both be uniax- 
ial dielectric-magnetic materials (with parallel symmetry axes) which are distributed as spherical particles 
[33]. Alternatively, four component materials could be homogenized: two isotropic dielectric materials and 
two isotropic magnetic materials, with all four component materials being composed of spheroidal parti- 
cles (all with the same alignment) [21]. However, the two-component formulation presented here is more 
computationally-efficient for the inverse Bruggeman formalism. 
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Figure 1: The constitutive parameters 711 (red, solid curve) and 722 (blue, dashed curve) plotted versus 
radial distance p. 
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Figure 2: A two-dimensional example of ray trajectories for pc = 1.1. Rays start at equally spaced 
locations along the line a;(0) = M ■ (30,1^,0) with —12.5 < v < 12.5 and where the rotation dyadic 
M = cos (/3 + yy) — sin ip {xy ~ y with 95 = 7r/4; and fc(0) is directed along (—1, —1, 0). 
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Figure 3: A three-dimensional example of ray trajectories for pc = 1.1. Rays start at equally spaced locations 
along the line x(0) = M • (30, v, 29) with —12.5 < v < 12.5 and where the rotation dyadic M is defined as 
in Fig. 2; and k{Q) is directed along (—1, —1, —1). 
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Figure 4: The constitutive parameter e;, = /i;, (red, solid curves), shape parameter [/f, = ?7a (green, dashed 
curves), and the volume faction /t, (blue, broken dashed curves) plotted versus radial distance p for (a) 
ta = \^a = 0.01, (b) Eq = = 0.1, and (c) da = Pa = 0.3. 
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